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ABSTRACT 

Consider  solutions  <H(x,e),G(x,e)>  of  the  von  Karman  equations  for  the 
swirling  flow  between  two  rotating  coaxial  disks 

1.1 )  eH^'^  +  HH'"  +  GG*  =0  , 

1.2)  eG"  +  HG'  -  H'G  =  0  . 

In  this  survey  we  describe  much  of  the  activity  of  the  past  30  years  - 
involving  physical  conjecture,  numerical  computation,  asymptotic  expansions 
and  rigorous  mathematical  results.  In  particular  we  focus  on  the  questions  of 
existence  and  nonuniqueness,  raonotonicity,  and  "scaling". 
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SIGNIFICANCE  AND  EXPLANATION 


Under  appropriate  conditions  the  steady-state  flow  of  fluid  between  two  planes 
rotating  about  a  common  axis  perpendicular  to  them  may  be  described  by  two  functions 
H(x,£),  G(X/e)  which  satisfy  the  coupled  system  of  ordinary  differential  equations 

+  HH"‘  +  GG'  =0 
EG"  +  HG'  -  H'G  =  0  . 

The  quantity  E  >  0  is  related  to  the  kinematic  viscosity  and  =  R  is  usually 
called  the  Reynolds  number. 

These  equations  have  received  quite  a  bit  of  attention.  First  of  all,  people 
who  are  truly  interested  in  the  phenomena  modeled  by  these  equations,  e.g.  fluid 
dynamicists,  are  Interested  in  this  problem.  However,  as  these  equations  have  been 
studied  by  a  variety  of  mathematical  methods,  they  have  taken  on  an  independent 
interest.  The  major  methods  employed  have  been  (i)  Matched  Asymptotic  Expansions 
and  (ii)  Numerical  Computations.  In  both  approaches  technical  problems  have 
appeared.  There  may  be  "turning  points,"  L.e.  points  at  which  H(x,E)  =  0,  Such 
points  require  special  and  delicate  analysis  within  the  theory  of  (i).  As  numerical 
problems,  these  equations  are  "stiff"  -  precisely  because  £  is  small.  The 
occurrence  of  "turning  points"  only  makes  computation  more  difficult. 

For  these  reasons,  these  equations  have  become  "test"  problems  for  methods  of 
"matching  in  the  presence  of  turning  points"  and  "stiff  O.D.E.  solvers."  However, 
when  one  has  "test  problems,"  one  needs  to  know  the  answers.  Unfortunately  here  the 
answers  are  largely  unknown. 

In  the  60  years  since  the  basic  von  Karman  paper  (1921)  and  the  30  years  since 
the  famous  papers  of  Batchelor  (1951)  and  Stewartson  (1953)  there  has  been  an 
intensive  interaction  between  physically  based  conjecture,  numerical  calculations, 
formal  asymptotic  expansions  and  rigorous  mathematical  results.  In  this  survey  we 
discuss  several  specific  questions  and  describe  this  interaction. 

This  article  is  an  expansion  of  a  talk  given  in  Oberwolfach,  Germany  at  the 
Conference  on  Singular  Perturbations. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  wit)i  the  author  of  this  report. 


incompressible,  steady  flow  -  "swirling  flow*.  Let  be  the  coordinates  of 

velocity  in  cylindrical  coordinates,  (r,d,x).  von  Karman  assumed  that  there  is  a  function 
H(x,e)  such  that 

q^  =  'H(x,£)  . 

Then  (see  (21,  [12])  there  is  a  function  G(x,e)  so  that  the  velocity  components  are 
described  by 

-  I  H*<x,e),  <J8  “  f  • 

The  functions  < H(x,e ),G( x,e ) )  satisfy  the  equations 

(1.1)  CK^''  +  HH'  "  +  GG‘  -  0  , 

(1.2)  EG"  +  KG'  -  H'G  -  0  . 

The  quantity  e  >  0  is  related  to  the  bulk  viscosity.  Equation  (1.1)  can  be  integrated  to 
yield 

(1.3)  eri"  •  ♦  HH"  +  G^  -  (H'  )^  =  U 

where  u  is  a  constant  of  integration. 

In  the  case  originally  studied  by  von  Karman,  the  flow  above  a  single  disk,  we  have  a 
problem  on  the  infinite  interval  (0,“)  and  the  constant  of  integration  is  known,  i.e.. 


where  »  G('“,e).  Moreover,  in  this  case  the  parameter  E  may  be  “scaled  out".  Assume 
*  0  and  let 

(1.4a)  e  »  x//e  , 

(1.4b)  H(x,C)  =  /r  h(C),  G(x,e)  =  g(e)G(0,E)  . 

Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG2q-80-C-0041 ,  and  by  the  Office 
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Then,  the  functions  <h((),g(C)  >  satisf/ 

<1.5)  h’*-  +  hh"  *  -j  ^  <h*)^  -  Y  ' 

(1.6)  q"  hg*  -  h’g  -  0  , 

and  the  boundary  conditions 

(1.7a)  h(0)  .>  0,  (no  penetration)  , 

(1.7b)  h'(0)  -  0,  (no  slip)  , 

(1.7c)  g<0)  -  1,  noraalizatlon  , 

( 1 . 7d )  h ( ( )  bounded  as  ^  , 

(1.7e)  g(C)  ♦  “  I*.'  t  ♦  “  . 

If  we  consider  the  flow  between  two  planes,  x  •  0,  x  >  1  rotating  with  constant 
angular  velocities  (*q/2»  (1j/2<  then  the  quantity  w  “  1i(e)  is  unknown.  This  latter  case 
was  first  studied  by  Batchelor  [2]  and  Stewartson  [38]  who  gave  conflicting  arguments  and 
conjectures.  In  this  case  the  boundary  conditions  are 
(I.Sa)  H(0,c)  e  H(1,e)  0,  [no  penetration] 

(l.ab)  H'(0,e)  »  H'()/£)  =  0,  Ino  slip) 

(1.8c)  G(O.e)  -  Oq,  G(1,e)  -  n^,  Ifljjl  +  |aj  *  0  . 

Both  of  these  problems  have  been  the  subject  of  many  numerical  studies  and  have  been 
attacked  by  formal  matched  asymptotic  expansions.  The  von  Karman  problem  was  studied 
numerically  by  0.  H.  Hannah  [8]  In  1947  and  by  H.  H.  Rogers  and  G.  N.  Lance  [35]  in  1962 
and  more  recently  by  D.  Dljkstra  and  P.  J.  Zandberger  [5]  in  1977  and  M.  Lentlni  and  H.  B. 
Keller  in  1980  (21).  These  recent  calculations  are  concerned  with  “traeing'out"  the 
branches  of  the  solution  set.  In  particular,  these  calculations  strongly  Imply  the  non¬ 
unicity  of  the  solution.  Formal  works  on  the  von  Karman  problem  have  been  carried  out  by 
M.  G.  Rogers  and  G.  N.  Lance  [35],  w.  G.  Cochran  (4),  H.  Ockenden  [28]  and  H.  K.  Kulken 
[191.  Rigorous  results  on  the  existence  and  uniqueness  question  are  not  complete.  J.  B. 
McLeod  (221  has  shown  existence  of  a  solution  for  all  non-negative  values  of  Q^.  In 
addition,  he  has  shown  non-existence  for  =  -1. 

For  the  two  disk  problem  numerical  calculations  have  been  carried  out  by  C.  E.  Pearson 
(30);  Lance  and  Rogers  (20);  I).  Greenspan  (7),  0.  Schultz  and  U.  Greenspan  [36li  G.  L. 
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Mellor,  P.  J.  Chappie  and  V.  K.  Stokes  (26);  N.  U.  Nguyen,  J.  P,  Ribault  and  P.  Plorent 
[271,-  S.  M.  Roberts  and  J.  S.  Shipman  (341;  H.  B.  Keller  and  R.  K-H.  Szeto  (13);  L.  O. 
Wilson  and  N.  L.  Schryer  (42);  G.  H.  Hoffman  (10);  H.  J.  Pesch  and  P.  Rentrop  (31);  M. 
Kubicek,  M.  Holodniok,  V.  Hlavacek  [11];  [17],  (18).  Formal  matched  asymptotic  expansion 
methods  have  been  used  by  A.  Watts  [40]  (who  also  did  numerical  calculations),  K.  K.  Tam 
(391,  H.  Rasmussen  [33],  B.  J.  Hatkowsky  and  W.  L,  Siegmann  [25].  Undoubtably  many  others 
have  also  worked  on  this  problem. 

As  in  the  case  of  the  single-disk  problem,  the  rigorous  mathematical  results  for  the 
two  disk  problem  are  incomplete.  The  basic  questions  of  ’existence*  and  ’uniqueness*  have 
remained  unanswered,  s.  P.  Hastings  [9]  and  A.  R.  Blcrat  [61  have  proven  existence  and 
uniqueness  for  large  £.  Their  arguments  are  essentially  a  perturbation  about  E  »  ».  ’J. 
B.  McLeod  and  S.  V,  Parter  [23]  considered  the  special  case  where  (1^  =  -0^  *  0.  They  have 
shown  the  existence  of  a  solution  for  all  E  >  0  and;  for  these  solutions,  they  gave  a 
complete  diaoussion  of  the  asymptotic  behavior.  More  recently  H.  O.  Kreiss  and  S.  V. 

Parter  [161  have  proven  the  existence  of  many  "large  amplitude"  solutions. 

Through  these  60  years  since  the  basic  von  Karman  paper  and  the  30  years  since  the 
Batchelor  paper  the  interaction  between  physically  based  conjecture,  numerical 
calculations,  formal  asymptotic  expansions  and  rigorous  mathematical  results  has  been 
intensive.  In  the  remainder  of  this  paper  we  will  discuss  several  specific  questions  and 
describe  this  interaction.  Of  course,  the  view  we  present  is  one  which  is  influenced  by 
our  own  work  and  interests. 

In  Section  2  we  discuss  the  counter-rotating  case:  flp  *  *  0.  In  Section  3  we 

discuss  the  monotone  co-rotating  case:  0  <  0^  <  ,  G'(x,E)  >  o.  Section  4  describes 

the  results  for  the  case  where  the  ’basic"  scaling  applies 
(1.9)  H(x,E)  =»  0(/e),  G(x,e)  .  0(1)  . 

Section  5  discusses  the  case  of  "order  i"  solutions.  In  Section  6  »*e  turn  to  the  question 
of  "cells".  Section  7  describes  the  existence  theory  for  "large  amplitude’  solutions. 
Finally  in  Section  8  we  discuss  some  more  unanswered  questions. 
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24  Counter-Rotating  Disks 

In  his  1951  paper  [2]  G#  K#  Batchelor  gave  special  attention  to  the  case 
(2.1)  G(0,e)  >  -1,  G(1,e)  .  1  . 

He  suggests  that  one  of  the  possible  solutions  the  main  body  of  the  fluid  would  be  In  two 
parts  with  different  angular  velocities  -  see  Figure  1  which  is  reproduced  from  [21.  In 
1952  K.  Stewartson  discussed  this  problem  and,  using  a  power  series  in  the  Reynolds  number 

R  =  1/6 

and  obtained  a  solution  in  which  the  core  has  (essentially)  zero  angular  velocity. 


O'letrtbubion  of  SLnqula.y  veloeitu  between  the  diece. 

Figure  1 

In  1965  C.  Pearson  (30)  computed  (numerically)  solutions  of  the  steady  state  problem 
as  the  (t  ♦  ■)  limit  of  a  transient  problem.  His  results  were  startling  in  that  his 
solutions  were  not  "odd"  about  x  «  V2  •  That  is 

H(x,e)  *  -H(1  -  x,e),  G(x,e)  *  ~G(1  -  x,e)  . 

Thus,  as  Serrin  (37)  observed,  Pearson's  results  implied  non-uniqueness.  Moreover,  the 
Pearson  solution  had  none  of  the  characteristics  suggested  by  either  Batchelor  or 
Stewartson  (see  Figure  2  -  taken  from  Pearson  (30)). 

In  1974  U.  H.  Hoffman  (10)  studied  this  problem  (among  otliers)  using  a  method  of 
computer  extension  of  the  Stewartson  perturbation  series. 
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Profiles  of  for  steady-state  motion  between  two  counter-rotating 

disks  with  R  ^  1000. 

Figure  2 

Tam  [391  (1969),  Rasmussen  [331  (1970),  and  Watts  [40)  (1974)  applied  matched 
asymptotic  expansions  to  this  problem. 

In  1974  J.  B.  McLeod  end  S.  V.  Parter  123}  proved  the  existence  of  an  "odd"  solution 


(odd  about 

X  * 

V2 ). 

Moreover,  they  gave 

a  complete  asymptotic  analysis  (e  '"0).  In 

particular. 

G(x 

,e) 

is  monotone,  i.e. 

(2.2) 

G’(x,e)  > 

0,  0  <  X  <  V2  • 

On  the  interval 

(  V2 

,1]  H(x,e)  is  characterized  by  three  points  x^,x2,X3  and  its 

negativity. 

We 

have 

(see  Figure  3) 

(2.3) 

H(x,e)  <  0, 

V2<  X  <  1  , 

(2.4a) 

H'(x,e)  <  0, 

V2  <  X  <  X,  , 

(2.40) 

H’ (x,e  )  >  0, 

X,  <  X  <  1  , 

-b- 


(2.5a) 

(2.5b) 

(2.6a) 

(2.6b) 


0  <  H*(x,e), 
H“(x,e)  <  0, 
0  <  H"'(x,£), 
H'  "(x,e)  <  0. 


<  x  <  Xj  , 
Xj  <  X  <  1  , 

<  X  <  Xj  , 


Xj  <  X  <  1  . 


Purtfiermore,  in  the  core  G(x,e)  is  exponentially  snail  while  in  the  >>oundary  layers  (at 
X  «  0  and  x  =  1)  the  solution  <  H(x,e),G(x,e) )  is  asymptotically  the  solution  of  a 
von  Karman  problem  with  =  0.  Finally,  consistent  with  the  remar)cs  above,  the  solution 
satisfies  the  basic  scaling  (1.9). 


Figure  3 


Thus,  123)  yielded  the  existence,  for  all  e  >  0,  of  an  odd  Stewartson  type 
solution.  Since  there  was  no  uniqueness  theorem,  these  results  did  not  exclude  the 
possibility  of  a  Batchelor  type  solution  or  a  Pearson  type  solution. 

One  result  of  [23]  asserts  that:  If  <  H(x,e ) ,G( x, e )  >  is  a  solution  which  is  odd 
about  X  =  V2  then  the  condition  (2.2)  implies  the  conditions  (2.3)-(2.7b)  and  vice- 
versa.  This  result  showed  that  the  computational  results  of  (7)  were  not  a  good 
approximation  to  a  solution  and  led  to  the  improved  numerical  method  described  in  [36).  In 
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addition,  the  singular  Batchelor  solution  Is  an  odd  solution.  His  picture  seems  to 


indicate  that  G'  >  0,  If  that  was  his  intention,  then  this  result  eliminates  the 
possibility  of  such  a  singular  solution.  Indeed,  since  Batchelor  expressed  his  equations 
in  such  a  way  as  to  imply  the  basic  scaling  (1.9),  the  singular  solution  is  ruled  out  by 
the  results  of  Section  4. 

Since  1974  there  has  been  further  computational  work  e.g.  [31 J  and  matched  asymptotic 
expansion  work,  e.g.  [251,  on  this  problem. 

Finally,  the  recent  results  of  Kreiss  anil  Barter  [16)  show  that  there  are  many  "large 
amplitude"  solutions  for  e  sufficiently  small.  It  seems  reasonable  that  the  Pearson 
solution  is  a  "single  hump"  large  amplitude  solution.  Watts  [401  came  to  this  conclusion 
on  the  basis  of  his  work  in  matched  asymptotic  expansions.  We  will  have  more  to  say  about 
this  in  Section  6. 
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3»  Co-Rotating  Disks 

Let  us  tlrst  consider  the  case 

(3.1)  ”o  ’  ”l  ”  ’’  i.e.,  GtO,e)  »  G(1,e)  -  l  . 

In  this  case  one  sees  at  least  one  solution  at  once 

(3.2)  H(x,e)  5  0,  G(x,e)  H  1  . 

It  is  not  difficult  to  see  that  this  solution  is  "stable",  i.e.  there  is  no  local 
bifurcation  -  see  [31. 

Let  G(1,c)  be  fixed,  i.e. 

(3.3)  G(1,e)  =  1  , 

and  let 

(3.4)  G(0,e)  =  s  . 

From  the  remark  above,  there  is  an  Interval  s^(£)  <  s  <  s^(E)  with  s^(e)  <  1  <  S2(e), 
so  that  there  is  a  solution  of  (1.1),  (1.2),  (1.8a),  (1.8b),  (3,3)  and  (3.4)  for 
3  6  (s^  (e  ),32(€ ) ).  This  remar)t  is  the  basis  for  a  numerical  method  studied  by  M.  Kubicek, 
M.  Holodniak  and  V.  HlavaSek  [18)  and  the  work  of  J.  Cerutti  [3).  However,  we  have  no 

knowledge  of  the  behavior  of  s^(e),  a^{c)  as  e  ♦  0+. 

The  formal  asymptotic  work  of  Stewartson  [38]  and  Watts  [40)  indicates  that  there  may 
be  other  solutions  of  this  problem  which  also  satisfy  the  basic  scaling  (1.9),  Watts 
suggests  there  is  a  solution  very  much  like  the  solution  obtained  for  the  counter-rotating 
case.  In  this  case  we  expect 

(3.5a)  H(x,e)  =  -H(1  -  x,e) 

(3.bb)  G(x,e)  =  G(l  -  x,e) 

with  G(x,e)  exponentially  small  in  the  core  and  H(x,e)  having  a  "shape"  similar  to  tne 
shape  described  by  ( 2. 3)-( 2. 6b) .  Some  of  our  preliminary  work  indicates  that  such 
solutions  do  indeed  exist  for  e  small  enough.  However,  we  conjecture  that  if  such 
solutions  exist,  then 

(3.6)  G(0,€)  <  0  . 
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In  t)ie  general  co-rotatinq  case 

(3.7)  0  <  G(0,e)  <  G(l,e)  . 

Batchelor  sui^yested  that  the  angular  velocity  would  he  monotone,  i.e. 

(3.8)  GMx,e)  >0,  0  <  X  <  1  . 

In  fact,  this  is  false.  McLeod  and  Parter  (241  have  shown: 

I^t  G(0,e),  G(l,e)  be  fixed  and  satisfy  (3.7).  Then  there  is  an  e  =  )  such  that 

if  0  <  G  <  £  and  if  there  is  a  solution  of  (t.1),  (1.2),  (1.8a),  (1.8b)  and  (i.Bc),  the 
ine<iuality  (i.B)  is  false. 

Once  more,  the  results  of  Kreiss  and  Parter  (16)  show  that  for  e  small  enough  there 
are  (many)  solutions  to  this  problem. 


4.  The  Basic  Scaling 

Many  of  the  authors  dealing  with  this  problem  have  assumed  the  basic  scaling  (1.9); 

(4.1)  |H(x,e)|  <  /e  B,  |G(x,e)|  <  B 

and  employed  the  change  of  variables  (1.4a),  (1.4b).  Thus,  one  has  the  equations  (1.5), 
(1.6)  on  the  larger  interval  [0,1//^1.  Prom  the  point-of-view  of  computation  this  leads 
to  regular  problems  -  albeit  on  a  large  interval  -  which  is  desirable  particularly  when 
employing  the  "shooting  method"  see  (201,  (261.  From  the  intuitive  point-of-view  as  well 
as  from  the  matched  asymptotic  viewpoint  it  is  reasonable  to  assume  that  the  solution  in 
the  boundary  layers  -  at  both  dislcs  -  behaves  li)ce  a  solution  of  the  von  Karman  problem 
which  then  "matches"  with  a  core  solution.  In  fact,  this  approach  was  used  by  Watts  (40], 
Tam  (391  and  Rasmussen  (331  -  and  implicitly  by  Stewartson. 

Within  this  context,  both  Batchelor  and  Stewartson  assumed  that:  in  the  core, 

6  <  X  <  1  -  S, 

(4.2)  q(5,£)  •  G(x,e)  ♦  G^,  a  constant  . 

They  considered  two  types  of  solutions 

Batchelor:  In  addition  to  (4.2)  we  have 

(4.3)  G.  *  0,  h>(5,e)  ♦  0  , 

i.e.,  the  core  rotates  as  a  rigid  body. 

Stewartson: 

(4.4)  G^  -  0  . 

Both  agreed  that  the  Batchelor  type  solution  would  appear  when 

0  <  G(0,c),  0  <  G(  1,E)  , 

i.e.,  the  co-rotating  case.  Stewartson  suggested  that  (4.4)  would  occur  when 

G(0,€)G(I,E>  <  0  . 

Solutions  of  this  type  have  been  obtained  both  numerically  and  via  matched  asymptotic 
expansions. 

In  fact,  if  the  basic  scaling  holds,  then  -  in  some  sense  (4.2)  must  hold.  This  fact 
IS  contained  in  the  following  results  of  Kreiss  and  Parter  (14). 
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Theorem  4.1  (see  Lemma  3.3,  Theorem  3.1  and  Theorem  4.1  of  [14]).  Let  4,  0  <  6  <  1/j  be 
given.  Let  (  H(  x,e^^ ) ,  G(x,e^ ) )  he  a  sequence  of  solutions  of  (1.1),  (1.2),  (1.8a),  (1.8b), 
(1.8c)  which  satisfy  (4.1)  for  some  constant  B.  Then  there  is  an  e(4)  and  an  M(4) 
such  that;  if  0<£  <e(6)  then  for  0<6<x<1-4<1  we  )iave 


1_ 

128 


(4.5)  I  ’AG  <x,e  )  -  u(e  )|  <  m(6)(i  +  B)e  "  *0  as  e  ♦  0+  . 

'  ^  n  n  '  n  n 

Obviously  one  can  extract  a  subsequence  e*  0+  so  that 

n 


(4.6a) 

(4.6b) 

Suppose  this  has  been  done.  If 
(4.7) 

then  there  is  a  constant  a  such  that 
(4.8a) 

In  fact,  both 


ti(e*  )  ♦  u  >  0  , 
n 


G(x,eM 

n 


♦  ±  2/? 


u  >  0  , 


H(x,e'  )//e^  ♦  a, 
n  n 


(4.8b) 


’/2G‘(x,e^)  -  ul. 


6  <  X  <  1  -  4 


-  H(x,E’) 

ru  " 

n 


are  exponentially  small  (in 

Remark;  Two  important  points  must  be  made.  First;  this  is  an  asymptotic  theorem,  there  is 
no  assertion  of  existence  of  solutions.  There  is  only  the  statement  that  if  such  solutions 
exist,  this  result  describes  their  asymptotic  behavior.  Secondly,  the  statement  that  in 
the  boundat y- layer  the  solution  is  essentially  the  solution  of  a  von  Karman  is  suggested 
but  this  discussion  is  not  entirely  complete. 

The  case  when  G^  =  0  is  more  complicated.  A  partial  discussion  is  given  in 
Section  6. 


I 
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5.  Order  One  Solutions 

While  the  basic  scaling  (1.9)  has  many  attractions  there  is  another  plausible 
scaling.  We  assume 

<S.t)  |h(x,£)|  +  |M’(x,e)|  +  lG(x,t)\  <  B  . 

If  this  bound  holds  then  the  physical  velocities  (q^.q^.q^)  are  bounded  in  any  cylinder 
r  <  R.  However,  in  order  to  guarantee  that  we  are  not  dealing  with  the  case  described 
earlier,  we  insist  that  H(x,c)  be  truly  of  order  1.  Specifically,  we  assume  there  is  a 
point  Xq,  0  <  Xg  <  1  and  a  constant  5  >  0  so  that 

(5.2)  |H(x^,e)|  >  6  >  0  . 

In  his  work  on  this  problem  with  matched  asymptotic  expansions  (33)  Rasmussen  had 
trouble  in  the  case  where  H  and  G  are  of  the  same  order.  It  has  been  suggested  that 
this  problem  is  involved  with  the  intrinsic  difficulties  of  Ac)cerberg-0'Malley  Resonance 
(11. 


In  fact,  the  matter  is  quite  simple:  (essentially)  there  are  no  such  solutions!! 
The  argument  is  in  two  parts.  If  (5.1)  holds  then  there  is  a  subsequence  ♦  0+ 
and  a  function  f)(x)  so  that 

(5.3)  H(x,e^)  ♦  8(x)  uniformly  on  (0,1)  . 

Further,  it  can  be  arranged  that  (5.2)  ta)jes  the  form 

(5.2 ' )  H(x  ,e' )  >  6  >  0  . 

V)  n 

Let  5  >  Xg  be  the  first  point  greater  than  Xg  at  which  H(B)  »  0.  Then 

(5.3)  H’<8)  -  0  . 


This  result  is  explicitly  given  as  Theorem  4.2  of  (15). 

It  now  follows  that  R(x)  >  0  and  is  of  the  following  form:  There  are 
®j>  ®  ^  <  •••  <  ”  1  and,  on  the  interval  '®j''’j+1*'  ^  " 

either  fl(x)  is  a  quadratic  or  H{x)  is  of  the  form 
(5.4a)  R(x)  »  A^(1  -  cosT^(x  -  0^)1  , 

where 


N  numbers 
0,1 ,. . .  ,N 


(5.4b) 


T , 
J 


2" 
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Finally,  while  this  result  is  never  explicitly  stateil  in  [15],  the  argument  given  in  [16] 


implies  that 
(5.5) 

Therefore,  (5. 


|G(0,eM|  +  lG(l,e|^)|  . 

)  is  a  necessary  condition  for  the  existence  of  "order  1"  solutions. 


I 

i 
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6.  Cells 


In  [26]  Hellor,  Chappie  and  Stokes  Introduced  the  concept  of  a  cell  and  computed 
several  multi-cell  solutions.  A  cell  is  the  region  between  successive  zeros  x^,X2  of 
H(x,e).  This  is  a  region  [in  (r,6,x)  space)  or  cell  In  which  a  portion  of  the  fluid  is 
"trapped*,  i.e.  the  fluid  cannot  cross  the  boundaries  x  -  x^ ,  x  -  X2.  Unfortunately  this 
definition  is  not  "tight*  enough.  It  allows  for  the  existence  of  cells  in  the  boundary 
layers  which  are  lost  as  c  *  0-t-. 

For  this  reason  we  have  adapted  the  following  approach:  Let  there  be  a  number  P  so 
that)  if 

(6.1a)  h(x,e)  =  e'^HCx.e), 

and  in  the  interior  of  (0.1) 

(6.1b)  h(x.e  )  ♦  R(x).  as  £  0  , 

n  n 

Definition:  A  "cell"  is  an  interval  (a. 8)  with  0  <  a  <  6  <  1  such  that 

(6.2a)  either  a  >  0  on  h(a)  <  0.  and 

(6.2b)  either  8  >  1  or  h(8)  ■  0.  and 

(6.2c)  |H(x)|  >0.  a<x<8. 

The  solution  obtained  in  (23)  in  the  counter  rotating  case  has  P  ■  -  and  leads  to 
two  oscillating  cells.  The  basic  result  is 

Theorem:  (see  Section  5  of  (15)).  Suppose  there  are  at  least  two  cells. 

.B,  )»  with 

(6.3a)  8^  < 

and.  these  cells  "oscillate*,  that  is 


(6.3b) 

fi(x)  >  0. 

0^  <  X 

<@1  . 

(6.3c) 

R(x)  <  0, 

^2  <  X 

<  02  • 

Then  R ( X ) 

has  most  4  cells.  Moreover. 

R(x)  is  a 

piecewise  quadratic  function  with  at 

most  two  distinct  pieces. 

That  is.  there 

is  a  point 

Xq  8  [0.1)  -  Note:  Xq  can  be  0 

or  1  -  and 

R(x)  has  the 

following  form 
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(6.4a) 


R(x)  e  c 


0  <  X  <  1 


(6.4b)  Ii(x)  «  +  ajX  +  83,  0  <  x  <  Xq  , 

(6.4c)  fi(x)  =  +  ^2*  *  **3'  ’'0^**1 

The  function  9(x)  la  exponentially  small  In  any  strict  Interior  subinterval  of  the 
two  Intervals  (O.Xq),  (xq,1). 

In  the  case  of  the  basic  scaling,  then  P  =  “  V2  •  Thus  we  see  that:  in  the  case  of 
-  Vo 

the  basic  scaling,  if  e  '2H(x,e)  is  convergent,  then  it  is  a  piecewise  quadratic  with  at 
most  two  pieces. 

Of  course,  Hellor,  Chappie  and  Stokes  used  the  basic  scaling  and  p  -  V2  i'*  their 

case. 

In  order  to  complete  the  discussion  of  the  basic  scaling  we  would  need  to  know  that 

-  Vo 

e  '2H(x,e) 

must  be  convergent  when  »  0. 
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7.  Large  Amplitude  Solutions 

In  140)  Watts  introduced  the  notion  of  a  large  amplitude  solution  < H(x, e ) ,G(x, e) >•  He 
obtained  a  formal  asymptotic  expansion  for  a  solution  with  one  "hump",  i.e.,  the  function 
H(x,e)  had  the  form 

H(x,e)  “  A(e)(1  -  co82iix1,  4  <  x  <  1  -  6  , 

where 

A(e)  ~  . 

The  function  G(x,e)  had  the  form 

G(x,e)  -  ±(2x)H(x,E)  . 

Watts  also  states  that  he  sees  no  reason  why  an  n-hump  solution  cannot  be  constructed  by 
the  same  method.  It  seems  he  also  obtained  the  relation  (7.6a)  (see  Figure  4  -  talcen  from 
Watts  (40)).  He  makes  no  mention  of  the  quantity  T. 

The  main  result  of  the  most  recent  paper  of  Kreiss  and  Parter  [16]  is  Theorem  7.1  (see 
Theorem  II  of  (16)).  Let  n  >  1  be  a  given  integer.  Let  s  be  a  fixed  real  number. 

Then,  there  is  an  e  •  e(n,8)  such  thatt  for  0  <  e  <  e  there  exists  a  solution 
<H(x,e),G(x,e)>  of  (1.1),  (1.2),  (1.8a),  (I.eb),  (1.8c)  with 
(7.1)  %  “  ’  • 

This  solution  may  be  described  as  follows.  There  are  exactly  (n  -f  1 )  numbers 

(7.2a)  0  «  a  (c)  <  0(e)  <  •••  <  o  (e)  <  o  (e)  -  l 

01  n-1  n 

at  which  H(x,e)  has  its  relative  minima,  i.e,. 


(7.2b) 

H'(o^(e),e)  »  0,  H"{o^(e),e)  >  0  . 

Between 

the 

0.  (c) 

J 

the  function  H(x,e)  is  essentially  positive.  That  is,  for 

4  >  0, 

24  < 

"j+l  ■ 

0.  we  have  -  for  small  e  - 
3 

(7.3) 

H(x,e)  >  0,  o^(e)  +  4  <  X  <  -  4  . 

The  numbers 

0  (c) 

] 

satisfy 

(7.4) 

where  T  is  a  fixed  number  which  will  be  described  in  the  Appendix.  The  number  is 

determined  from  this  relationship  and  (7.2a).  We  have 
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The  function  G(x,e)  has  at  least  n  nodal  zeros:  0  <  Y. (£)  <  Y,(e)  <  •••  <  Y  (t)  <  !• 

t  Z  n 


If  s  ^  0  and 


Yj(e)  <  <i^(e),  a^(e)  -  r^(e)  »  0(e)  . 


sign  s  =  (-1)*^ 


Then  G(x,e)  has  (n  +  1 )  zeros.  The  additional  zero,  Tq(6)  satisfies 


Finally 


0  <  Ygte)  -  0(E)  . 


IHI  -  e~^,  IGI  -  e'^  , 


and,  on  the  interval  a.  +S<x<0.^, -6  we  have 

]  3+1 

(7.7b)  e^G(x,e)  -  e^(-i )^*’t^H(x,e )  , 

(7.7c)  e^H(x,e)  "  Aj(l  -  cose^(x  -  o^))  . 

Before  going  on  to  a  sketch  of  the  proof,  it  is  worthwhile  to  see  what  are  some  of  the 
consequences  of  this  theorem,  when  discussing  "branches*  of  solution  of  homogeneous  second 
order  equations  it  is  useful  to  characterize  solutions  by  the  number  of  interior  zeros. 
Since  a(x,e)  satisfies  such  a  homogeneous  second  order  equation,  let  us  characterize  the 
solution  pair  <  H(x,c ),G(x,e ) )  by  the  zeros  of  G(x,e), 

Case  1 :  s  >  0.  For  every  even  n  >  2  there  are  at  least  two  solutions  <  H,G>,  <  H,G> 
with  G(x,e ),a(x,c )  having  exactly  n  interior  zeros  and 


H(x,e)  >  0,  H(x,e)  >  0,  (essentially)  . 


Proof !  I,et  n  =  n.  Since 

sign  8  -  (-1 )"  >  0 
sign  3  *  (-l)"*^ 

the  function  G(x,e)  of  the  solution  pair  <H,G>  described  in  Theorem  7.1  has  exactly 
n  interior  zeros.  However,  the  function  G(x,e)  of  the  solution  pair  (  H,G)  associated 


with  n  -  n  -  1  also  has  exactly  n  •  (n  -  1 )  +  1  interior  zeros. 

Case  2:  s  <  0.  For  every  odd  n  >  3  there  are  at  least  two  solutions  (  H,G>,  (  H,G  >  with 


G(x,e),  G(x,f. )  having  exactly  n  interior  zeros  and 


H(x,e)  >  0,  H(x,c)  >  0,  (essentially)  . 
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If  n  =  1  there  is  at  least  one  solution  ( H,G >  with  G(x,e)  having  exactly  one  interior 
zero  while  H(x,e)  is  essentially  positive. 

Proof:  Let  n  =  n.  Then 

-1  “  sign  s  *  (-1)"'*'’  •  1  . 

Then,  for  the  solution  < H,G  >  described  in  Theorem  7.1  G(x,e)  has  exactly  n  =  n 
interior  zeros.  If  n  >  3  then  let  n  =  n  -  1.  Since 

-1  =  sign  s  =  (-1)"^’  =  (-1)" 

the  function  G(x,e)  of  the  solution  H,G  associated  with  n  ^  n  -  1  has  exactly 
n  +  1  =  n  interior  zeros. 

Case  3:  s  =  0.  For  every  n  >  1  (even  or  odd)  there  is  at  least  one  solution 
(  H(x,e ),G(x,e ) >  with  G(x,e)  having  exactly  n  interior  zeros  while  H(x,e)  is 
essentially  positive. 

Proof i  Let  n  =  n  and  let  <  H,G >  be  the  solution  described  in  Theorem  7,1. 

This  theorem  is  proven  by  first  obtaining  an  0(1)  solution  with  boundary  values 
G(0,e),  G(1,€)  which  are  0(€^'^^).  These  "pathological"  solutions  are  obtained  from  a 
"shooting"  argument. 

The  arguments  and  results  of  (IS)  are  used  to  show  that  when  H(x,e)  >  then 

( H(x,e ),G(x,e ))  must  have  the  form 


(7.8a) 


H(x,e) 


cos  T ( X  -  o ) )  , 


(7.8b)  G(x,e)  "  TH(x,e)  . 

On  the  other  hand,  when 

(7.9)  H(x,e)  .  0(e^/^) 

one  employs  the  change  of  variables 


(7.10a) 


5 


^  -  *0 
t’/3 


(7.10b) 


h(e,e)  =  e"^''^H(x,e),  g(«,c)  -  e"^^^G(x,e)  . 
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The  functions  (  h,g  >  satisfy  the  etjuations 

(7.na)  n'"  +  hh"  +  ^  -  1  (h‘)^  -  , 

(7,110)  g*  +  hg'  -  h*g  »  0  . 

The  initial  values  are  chosen  so  that 
(7.12)  |li|  <  kc,  lu/e^^^l  <  . 

Now  it  is  not  difficult  to  see  that  h(C,k)  converges  to  a  quadratic  function  h(6)  of 
the  form 

(7.13a)  h(5)  « -j  hji?  -  )^«  hj  a  constant  . 

Thus.  g(C,t)  converges  to  a  function  g(C)  which  satisfies 
(7.13b)  g"  +  Hg*  -  H’g  -  0  . 

The  final  result  depends  on  an  elementary  degree  theory  argument  and  an  analysis  of 
the  solutions  of 

(7.14)  g-  +  e^i*  -  2Ci  -  0  . 

The  facts  about  this  equation  are  described  in  the  Appendix. 


8.  Coamants  and  Questions 

Despite  all  we  now  know  about  the  solutions  of  (1.n>  (1.2),  (1.8a),  (1.8b),  (1.8c), 
there  are  many  interesting  questions  still  unanswered. 

Question  1 ;  It  is  not  difficult  to  see  that  the  solutions  obtained  by  Hastings  (91  and 
Elcrat  (6]  for  large  e  >>  1  can  be  obtained  by  an  iterative  procedure  with  e  as  a 
parameter.  Thus,  for  e  large  we  have  a  curve  of  solutions. 

The  first  question  is:  keeping  fixed  and  varying  e  downward,  how  far  can 

these  continue  of  solutions  be  continued?? 

It  is  reasonable  to  assume  that  these  solutions  exist  for  all  e  >  0.  The  next 
question  is;  if  that  is  so,  do  these  solutions  < H,G >  satisfy  the  basic  scaling  (1.9), 
i.e. 

(8.1  )  H  »  0(/e),  G  =  0(1  )  ?? 

Question  2;  With  fixed,  are  there  families  of  solutions  which  satisfy  (8.1)?? 

Recall  that  when  (1^  -  -0^  the  solutions  obtained  by  McLeod-Parter  (231  do  indeed  satisfy 
(8.1).  And  of  course,  the  trivial  solution  (3.2)  satisfies  (8.1).  The  computational 
evidence  suggests  that  there  are  such  solutions. 

Question  3;  If  there  are  solutions  satisfying  (8.1),  are  they  unique??  The  computational 
evidence  suggests  that  there  are  many  solutions.  If  there  are  many  solutions  -  how  does 
one  characterize  the  possible  values  G^  of  (4.2).  The  results  of  Rasmussen  (331  and  the 
conjecture  of  Stewartson  ( 381  suggests  that  whenever 

0,  <  0 

we  must  have 

G^  »  0  . 

Question  4;  Do  the  large  amplitude  solutions  lie  on  continua  of  solutions?  In  particular, 
given  s  and  n  >  2,  is  there  a  continuum  of  solutions  which  exist  for  e  <  G(n,8)  which 
contains  the  pair  <H,G>,  <H,G)  of  solutions  with  G(x,e)  and  G(x,e)  having  exactly  n 
interior  zeros  described  in  Section  ??? 

giestion  5;  Are  there  "other"  solutions??  That  is;  are  there  families  of  solutions  which 

—2  »2 

do  not  satisfy  (8.1)  other  than  those  large  amplitude  solutions  (H  ”  e  ,G  ~  e  )  found 
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in  (16)??  Of  course,  for  each  solution  -H(1  -  x,e),G(1  -  x,e)  Is  again  a  solution. 
Hence  there  are  related  large  amplitude  solutions  with  H(x,e)  <  0  (essentially). 

Question  6:  The  function  (first  discussed  by  McLeod  [22]) 

♦(x.e)  -  t(G')^  +  (H“)2) 

plays  an  important  role  in  much  of  the  analysis  -  see  [22],  [23],  [14].  It  is 
characterized  by  the  fact  that  there  Is  a  unique  point  Y  ~  Y(c)  6  [0,1]  at  which  ^(x,e) 
has  a  minimum.  Furthermore,  if  4(x,e)  has  a  relative  minimum,  it  is  the  minimum.  Watts 
[40]  observes  that  -  asymptotically  -  for  the  'solutions*  which  he  obtain  which  also 
satisfy  (8.1), 

(8.2)  0  <  Urn  Y(e)  <  1 
while  in  the  case  of  his  large  amplitude  solutions 

(8.3)  Lim  Y(e)  -  0  or  1  . 

It  is  not  difficult  to  show  that  the  large  amplitude  solutions  constructed  in  [16)  do 
indeed  satisfy  (8.3).  The  question  is:  does  Lim  Y(t)  characterize  the  "size”  of  all 
solutions? 

Question  7;  The  negative  results  of  McLeod  and  Farter  in  Section  3  and  the  monotonicity  of 
G(x,e)  for  the  counter-rotating  case  leads  to  the  following  observation  and  question. 

For  flg  ■  «  -1  and  all  e  >  0  there  is  a  solution  <H,G>  with 

G'(x,e)  >  0  . 

For,  0  ^  that  statement  cannot  be  true.  Therefore,  the  question  Is:  if  there  is 

a  number  a,  -1  <  a  <  0  determined  such  that;  if  »  1 ,  and 
(8.2)  -1  <  <  a 

then,  for  all  e  >  o  there  is  a  solution  <H,G)  with  G'(x,e)  >  o.  if  there  is  no  such 
number,  then  -  -1  is  a  very  special  case  Indeed. 

The  final  question  is  a  very  large  one.  Given  a  solution  how  can  one  determine  its 
time-dependent  stability? 
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Appendix 

A  key  part  of  the  argument  leading  to  results  of  Section  7  -  l«e«  the  results  of  [16] 

•  is  the  analysis  of  the  solutions  of 

2 

(A.l)  g"  +  6x  g*  -  26xg  »  0  • 

We  sketch  this  analysis. 

We  can  restrict  ourselves  to  the  case  6  »  1.  Let  g(x;1)  be  a  solution  of  (A.I) 
with  6=1,  Then,  for  any  6  >  0,  a  direct  calculation  shows  that 
(A. 2)  Y(x;6)  =  g(S'^^x;1) 

Is  a  solution  of  (A. 1)  with  this  value  of  6. 

For  the  remainder  of  this  discussion  we  have  6  -  i . 

Using  the  W.K.B.J.  method  (see  Chapter  6  of  [29]  and  the  method  described  by  Wasow  In 
[41,  pp.  52-61]  we  see  that  there  are  two  linearly  Independent  solutions  g^(x),  gjtx)  sod 


(A.  3a) 

91 (x)  ^ 

x^(l  + 

2/3x2), 

X 

(A. 3b) 

q;(x)  - 

2x(1  ♦ 

2/3x2), 

X 

(A. 3c) 

9?  -  2, 

X 

3 

(A. 3d) 

92(x)  ' 

X  ^expl 

X 

Similarly,  there  are  two  linearly  Independent  solutions  <^^(x), 


(A. 4a)  V’i(x)  '  X  ^exp[- x  +  +"  , 

(A. 4b)  X  ■>•  +“  . 

Since  the  function  g)(x)  can  be  written  as  a  linear  combination  of  i^i(x)  and 

see  that  there  is  a  unique  constant  T  such  that 

(A. 5)  g,(x)/x2  ♦  T,  X  +  +-  . 

Of  course,  this  x  Is  the  quantity  of  Section  7. 

We  don't  need  to  hnow  much  about  the  functions  i^^(x),  <^2(ix).  It  suffices  that 
(A. 6)  V’jCO)  *  0  . 

This  elementary  result  follows  almost  immediately  from  the  maximum  principle  or  the 
representations 
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(A.  7a) 


{'<’,exp(^]|  =  2x'<>^exp[-^]  , 

,  .  3  3 

(A. 7b)  ^  {'0,exp[^l}  =  2v!^exp[^]  , 

and  the  fact 

(A. 8)  Vfljco)  =  0  . 

Our  major  Interest  centers  on  the  function  g^(x).  The  basic  facts  are: 

(A. 9a)  g'(x)  <  0,  -•  <  x  <  •  . 

There  is  a  value,  say  g,  at  which  g^(§)  =  0-  This  unique  zero  can  be  estimated  by 
(A. 9b)  -1  <  g  <  0  . 

Finally 

(A. 9c)  T  <  0  . 

These  results  are  obtained  by  a  detailed  argument  based  on  elementary  considerations, 
the  maximum  principle,  the  oscillation  theorem  and  the  series  expansion  of  the  two 
functions  Y^(x),  which  satisfy  (A.1)  -  with  5=1-  and  also  satisfy 

(A. 10a)  Y,(0)  =  0,  yJ(0)  -  -1  , 

(A. 10b)  “  0  • 

We  remar)t  that  it  is  equally  easy  to  obtain  these  results  by  rigorous,  careful, 
numerical  computation.  In  fact,  computations  by  Jerry  Browning  of  NCAR  indicate 

T  =  -2  . 
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